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ABSTRACT 



Context: The expansion of network magnetic fields with height is a fundamental property of flux tube models. A rapid expansion is 
required to form a magnetic canopy. 

Aims: We characterize the observed expansion properties of magnetic network elements and compare them with the thin flux tube and 
sheet approximations, as well as with magnetoconvection simulations. 

Methods: We used data from the Hinode SOT NFI NaDi channel and spectropolarimeter to study the appearance of magnetic flux 
concentrations seen in circular polarization as a function of position on the solar disk. We compared the observations with synthetic 
observables from models based on the thin flux tube approximation and magnetoconvection simulations with two diff'erent upper 
boundary conditions for the magnetic field (potential and vertical). 

Results: The observed circular polarization signal of magnetic flux concentrations changes from unipolar at disk center to bipolar near 
the limb, which implies an expanding magnetic field. The observed expansion agrees with expansion properties derived from the thin 
flux sheet and tube approximations. Magnetoconvection simulations with a potential field as the upper boundary condition for the 
magnetic field also produce bipolar features near the limb while a simulation with a vertical field boundary condition does not. 
Conclusions: The near-limb apparent bipolar magnetic features seen in high-resolution Hinode observations can be interpreted using 
a simple flux sheet or tube model. This lends further support to the idea that magnetic features with vastly varying sizes have similar 
relative expansion rates. The numerical simulations presented here are less useful in interpreting the expansion since the diagnostics 
we are interested in are strongly influenced by the choice of the upper boundary condition for the magnetic field in the purely 
photospheric simulations. 

Key words. Sun: magnetic topology - Sun: photosphere 



1. Introduction 

At present, observations do not provide an explicit picture 
of how the chromospheric network magnetic field is struc- 
tured. On one hand, we have increasing observational evi- 
dence of something which can be interpreted loosely as a 
canopy structure: e .g., fibrils in the Ca II infrared triplet lines; 
IVecchio et al.ll2007i large-scale canopy structures in combined 
Zeeman and Hanle studies: |Biandaetal. (1998); Stenfloetal. 
(12002), canopy-like expansion seen in magnetograms near the 
limb: lJones & Giovanellil (11983.) ). On the other hand, we know 
that the Sun is more complicated than implied by simple mod- 
els (IWed emeyer- Bohm et al.ll20()9l) . An important question is to 
what extent simple models of flux tubes are able to reproduce the 
center-to-limb appearance of the network magnetic field struc- 
tures. 

The appearance of network magnetic flux concentrations in 
circular polarization maps changes from unipolar at disk center 
to bipolar near the limb. This is consistent with network mag- 
netic fields expanding and fanning out with height as proposed 
by iGabriell (Il976|). In the chromosphere , the fanning is clearly 
present; e.g.. I Jones & Giovanellil (Il983l) found low-lying, 200- 
800 km, magnetic can opies in magnetog rams taken near the so- 
lar limb. More recently 'Kontar et aD (|2008) has used hard X-ray 
observations from RHESSI to estimate the expansion and found 
that the magnetic field expanded noticeably at a height of ^ 900 
km. 



Expansion of magnetic field with height has been studied in 
photospheric structures mostly using magnetic flux tube models. 
Indeed, flux tube models predict a rapid expansion of the field 
with height. Solanki et al. (1999) shows that magnetic structures 
as diflTerent in size and flux as small flux tubes and sunspots have 
similar relative expansion rates, which agree with the thin flux 
tube approximation. A study of the characteristics of magnetic 
flux structures in radiative magnetohydrodynamic (MHD) simu- 
lations revealed the expansion properties to be similar with the 
thin flux tube and sheet approximations ( Yelles Chaouche et al.l 
2009). The expansion is seen in observations: e.g., a thin flux 
tube model can simultaneously reproduce the observed Zeeman 
splittings of Mg 12.32 yum, Fe 525.0 nm and Fe 1.56 yum lines, 
which span the upper to the lower photosphere in formation 
height (lBruls&Solankilll995l) . Bruls & Solanki also showed 
that a flux tube model can explain the Mg 12.32 yum line pro- 
file shapes observed by Zirin & Popp (1989). Additional evi- 
dence of expansion is that a canopy resulting from the expan- 
sion of a flux tube can best explain the observed photospheric 
asymmetric Stokes V profiles with weak zero-crossing shifts 
(iGrossmann Doerth et aDll988l) . 

In this paper we use circula r polarization maps from the 
Solar Optical Telescope (SOT, iTsuneta et aP l2008h on the 
Hinode satellite to study the average expansion properties by 
characterizing how the appearance of network flux concentra- 
tions changes from the solar disk center to the limb. The center- 
to-limb approach lets us examine the expansion at diflferent 
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viewing angles and at different heights due to the shift in the 
height range where spectral lines are formed as a function of [i 
(jd = cos(0), where 6 is the viewing angle). To further expand the 
coverage we use SOT observations from the spectropolarimeter 
(SP) and the narrowband filter imager (NFI) NaDi channel. We 
combine the observations with modeling the expansion of mag- 
netic flux with height at various yu-values by using the thin flux 
tube and sheet approximations and more realistic 3 -dimensional 
magneto-convection simulations. 



2. Observations and models 

2.1. Observations and analysis 

The SOT NFI NaDi (589.6 nm, efl'ective g factor 1.33) data 
used in this study consist of 19 circular polarization filtergrams 
of non-active region maps at diflTerent positions on the solar 
disk (see Table [T] for time, location, exposure time and field 
of view, where xcen and ycen are orthogonal distances from 
Sun center and xcen is directed along the equator). The data 
reduction was done using the Solar Software (SSW) Q pack- 
age routine fg-prep.pro. To increase the signal-to-noise ratio, 
frames taken within 3 minutes are summed together (number of 
frames is shown as a multiplication factor in from of the ex- 
posure time in Table [T]). Errors in pointing are corrected in the 
limb data sets by forcing yu=0 to coincide with the visible so- 
lar limb. The NaDi line is often called a chromospheric diag- 
nostic but 3 -dimensional non-local thermodynamic equilibrium 
(nonLTE) radiative transfer calculations of magnetoconvection 
simulations have demonstrated that the line is to a large de- 
gree photospheric, es pecially in magnetic flux concentrations 
(iLeenaarts et al.ll2QQ9h . 

For SP data we use 10 raster scans with exposure times of 
eight seconds or more (i.e., deep magnetogram mode) at varying 
positions on the solar disk (see Table [B reduced using the Solar 
Software package routines sp-prep.pro and stksimagessbsp.pro. 
The step size is 0.149 '' in all scans but number 5 where it is 
0.29 '' . The reduced data consist of maps of circular polarization 
and linear polarization in the 630.1 (eflTective g=l.61) and 630.2 
nm (eflTective g=2.5) Fe I l ines, both combine d and separate. For 
details of the routines see iLites et al.l (l2008l) . As in NaDi data, 
pointing errors are corrected in maps where the solar limb is 
visible. 

Flux concentrations (coherent features with circular polar- 
ization signal clearly above the noise level) are identified manu- 
ally in each NFI and SP circular polarization image. For the iden- 
tification we use images of the absolute value of the flux instead 
of the signed flux to reduce the bias towards choosing bipolar 
features near the limb. A central pixel of each identified feature 
is chosen and a vector connecting the disk center and limb pass- 
ing through this point is defined. This vector is the average radial 
cut of the feature. The SP circular and linear polarization radial 
cuts are made by averaging over 7 pixels perpendicular to the 
vector for the 630.1 and 630.2 nm. Unless otherwise specified 
we use the average of the two lines. For NFI data the radial cuts 
(width 7 pixels) are made of the Stokes V filter signal. No cor- 
rection for foreshortening is made. Finally, peaks (location and 
amplitude) in the radial cuts are identified manually in each cut 
(see Fig.[T]). Since we are interested in the large scale expansion 
of the field as opposed to a detailed study on the small scale mor- 
phology of the structures, spatial and, in the case of NaDi also 
temporal averaging, are justified. 



http://www.lmsal.com/solarsoft 



2.2. Models 

We have employed two independent sets of models, the sim- 
ple thin-tube approximation and realistic 3-D simulations. The 
thin flux tube model is based on the assumption that the flux 
tube is much thinner than the pressure and density scale heights. 
The plasma and magnetic properties inside the t ube can then be 
expanded in terms of a low order polynomial (lDefouwlll976l: 
[Roberts & Webb|[l978h . Here we use the lowest order polyno- 
mial possible: 0^^ order in radius, r, for the vertical component 
of the magnetic field, B^, and V^ order in the radial component of 
the field, Br. The component of the Lorentz force corresponding 
to magnetic tension does not appear in this 0^^ order in r expan- 
sion of Bz (it first appears in the 2"^ order expansion). We use 
semi-empirical atmospheric models for the internal and external 
models, from which the magnetic field can be derived using the 
condition that the total (gas -h magnetic) pressure inside the tube 
at each height is equal to that of the external gas pressure. 

The external non-magnetic atmosphere is given by the semi- 
empirical SRPM mo del of the quiet Sun low chromosphere 
(i Fontenla et al .1120071) and internal atmosp here by the recent fac- 
ular and plage models of the photosphere dFontenla et al.ll2006l) 
which have been extended in height by adding the FALF (facular 
model. iFontenla et aPl 19931) model chromosphere. The eff'ect of 
internetwork fields is not considered in the external model. This 
is justified since their strength in the lower photosphere is on the 
order of a 100 G, i.e., the magnetic pressure is roughly 1% of the 
gas pres sure. At greater heights the se internetwork fields fall off* 
rapidly (ISchtissler & Vogle^l2008h and so remain unimportant 
to this particular study. We produce 2-dimensional atmospheric 
slabs (x,z) by varying the flux tube/sheet parameters: radius at 
z=0 km (ro), field strength at z=0 km (^o), and the internal flux 
tube/sheet atmosphere (facular or plage). 

Observables are synthesized under t he assumption of LTE 
using the Nicole-code of Socas-Navarro (!Socas-NavarrQll2QQl1) . 
Arguably, especially the NaDi line is influenced by nonLTE ef- 
fects, but since we consider the Stokes V signal which has its' 
maximum value away from the line core (where nonLTE ef- 
fects are significant) and , more specifically, in the ratio of the 
Stokes V amplitudes, we chose to do the synthesis in LTE. We 
are not interested in comparing absolute values of the ampli- 
tudes or details of the line profiles, which may well be aflfected 
by nonLTE eff'ects, but rather in the ratio between the ampli- 
tudes of the center- and limb- sides of magnetic concentrations, 
— , which should be independent of nonLTE eflfects, at least to 
first order. We compute 1 -dimensional rays passing through the 
slabs at various line of sight (LOS) angles and produce synthetic 
observables. Synthetic cuts are created from the Stokes V ampli- 
tudes and areas (sign determined by the sign of the blue lobe). 
This is done also for the NaDi line instead of using the filter 
transmission profile. By doing this the NaDi synthetic observ- 
ables pick up the maximum signal at each point. We compute 
the observables, radial cuts of Stokes V areas and amplitudes, 
using the original tube/sheet model resolution (10 km in the hor- 
izontal direction, rays placed at 10 km intervals) and also for the 
spectral (90 mA for NFI, 23 mA for SP) and spatial (theoreti- 
cal point spread function (PSF) for a 50 cm telescope, 0.08'' and 
0.16'' pixels size for NFI and SP, respectively) resolutions com- 
parable to those of Hinode. For the purposes of this study, we 
ignore contributions from the spider and other optical elements 
of Hinode on the PSF. Additionally, we assume the PSF to be 
1 -dimensional. These idealizations are appropriate since we are 
not interested in the detailed structure of the features in our sim- 
plified model. 
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Table 1. NFI and SP data sets 



NFI 
















date 


xcen [''] 


ycen ["] 


^ 


exp [s] 


FOV V X "] 


1 


2008-01-16X11:13:34.821 


4 





1.00 


4 X 0.20 


113. X 113. 


2 


2008-04-03X12:55:37.304 


100 


-270 


0.95 


4 X 0.20 


64.0 X 164. 


3 


2007-11-11X13:07:55.041 


-313 


-57 


0.94 


10 X 0.20 


41.0 X 164. 


4 


2007-12-01X11:04:34.300 


325 


3 


0.94 


4x0.12 


225.x 113. 


5 


2007-11-29X11:28:37.095 


377 


-184 


0.90 


4 X 0.20 


113.x 113. 


6 


2007-11-05X05:28:45.378 


503 


109 


0.84 


10 X 0.20 


41.0 X 164. 


7 


2007-11-05X07:27:33.212 


518 


110 


0.83 


10 X 0.20 


41.0 X 164. 


8 


2008-05-13X16:27:02.821 


511 


187 


0.82 


7 X 0.20 


56.3x113. 


9 


2007-11-15X10:09:35.952 


-58 


-561 


0.81 


10 X 0.20 


41.0 X 164. 


10 


2007-11-15X08:23:31.692 


-71 


-561 


0.81 


10 X 0.20 


41.0 X 164. 


11 


2008-04-05X20:34:34.811 


542 


-291 


0.77 


4 X 0.20 


64.0 X 164. 


12 


2007-11-13X00:44:33.710 


-303 


-551 


0.76 


10 X 0.20 


41.0 X 164. 


13 


2008-05-14X12:51:32.612 


735 


-60 


0.64 


4 X 0.20 


64.0 X 164. 


14 


2007-11-17X10:45:35.227 


858 


-114 


0.43 


3 X 0.20 


164. X 164. 


15 


2008-03-24X14:28:38.464 


-10 


-947 


0.31 


4 X 0.20 


64.0 X 164. 


16 


2008-03-22X16:48:40.335 


-10 


-946 


0.31 


4 X 0.20 


64.0 X 164. 


17 


2007-11-10X11:06:02.179 


-6 


928 


0.25 


7 X 0.20 


164. X 164. 


18 


2007-11-12X11:37:31.448 


230 


912 


0.20 


7 X 0.20 


164. X 164. 


19 


2007-12-19X09:34:40.109 


931 


-163 


0.19 


4x0.15 


113.x 113. 


SP 


1 


2007-09-06X15:29:04.938 


-46 


6 


1.00 


8.0 


149 X 162 


2 


2007-09-10X15:03:12.528 


87 


7 


1.00 


9.6 


56 X 162 


3 


2007-09-22X09:28:05.954 


27 


307 


0.95 


13. 


47 X 162 


4 


2007-09-21X11:28:34.622 


25 


607 


0.77 


13. 


47 X 162 


5 


2007-09-10X22:13:23.366 


-247 


656 


0.68 


1.6 


153 X 162 


6 


2008-01-21X06:02:48.747 


760 


-387 


0.46 


13. 


56 X 162 


7 


2007-09-07X03:59:46.805 


-18 


922 


0.28 


9.6 


105 X 162 


8 


2007-09-03X16:06:00.480 


144 


922 


0.23 


9.6 


298 X 162 


9 


2007-11-03X05:17:25.426 


54 


952 


0.11 


13. 


208 X 162 


10 


2007-11-01X06:18:20.542 


93 


952 


0.080 


13. 


248 X 162 



For the purposes of this paper the difference between flux 
tubes and sheets lies in their linear expansion rates. In the thin 
tube approximation the expansion is driven by the magnetic flux 
conservation with height: BA = const, where B is the field 
strength and A is the tube's cross-section. As B drops roughly 
exponentially with height, A increases inversely. In a flux tube 
the increase in radius goes as r ~ Va, while in a sheet the 
magnetic field only expands in the direction perpendicular to 
the sheet, so that r ^ A. Consequently, r increases far more 
rapidly in sheets than tubes. The ratio, ^, strongly depends on 
the viewing angle. It is smallest when rays are perpendicular to 
the sheet direction and infinite, i.e., unipolar, when rays are par- 
allel. Therefore, to delimit the range we only synthesize radial 
cuts for rays pointing in the perpendicular direction. The expan- 
sion rate for a magnetic feature consisting of several flux tubes, 
which merge above a height determined by the density of the 
features and their magnetic flux follows the same scaling. Below 
the merging height the flux tubes expand individually. Above 
this, they expand as a single feature following A(z) ~ Ao/B(z), 
where Aq is the summed cross-section of the individual tubes at 
z = 0. 

To study flux expansion under more realistic conditions 
we use radiative magnetohydrodynamic (MHD) simulations. 
In contrast to the simple flux tube model, the simulations al- 
low us to look in detail at the eflfect of small scale, inho- 
mogeneous structures on the expansion. We use the MuRAM 
code (MPS/University of Chicago Radiative MHD, Vogler et al. 
l2005h . The code takes into account eff'ects of full compressibil- 
ity, open boundary conditions as well as non-gray radiative trans- 
fer and partial ionization. In the current study the simulation do- 



main is 24Mmx24Mmxl.68Mm (576x576x120 grid points) of 
which approximately 700 km is above the logirsoonm) = level. 
To construct our initial condition we began with a snapshot from 
a non-magnetic convection run. To this we added a magnetic 
strip with a squared Gaussian profile. We chose a peak value 
of 300 G so that the initial magnetic field is energetically im- 
portant above the surface but much less so below the surface. 
The full width at half maximum was chosen to be 3.157 Mm, 
so that the patch of magnetic field spans several granules. This 
setup is reminiscent of an isolated, unipolar network lane. The 
field then evolves on several timescales. The slowest of these is 
the dispersion of magnetic flux due to the granular motions ad- 
vecting them. This process can be modeled as a random walk 
with a length- scale of approximately 1 Mm and a timescale of 
approximately 5-10 minutes. A second important timescale is 
that required to achieve pressure balance in the vertical direc- 
tion and for our box, which extends 700 km above approxi- 
mately log(T5oonm) = surfacc and where the sound-speed is 
^ 10 km/s, is some small multiple of 700 km/Cs=70 s. The time- 
scale for pressure balance to be achieved in the horizontal direc- 
tion should be similar because for fields of substantial strength, 
the length scales of the fields in the horizontal and vertical direc- 
tions are likely to be similar. An intermediate timescale is that 
of the magnetic flux to be expelled from the interior of the gran- 
ules, which takes place on the convective overturning timescale 
of the flows. Since we study the expansion with height of a thin 
strip of magnetic field, we allowed the simulation to evolve from 
the initial condition for a few granular lifetimes (^15 minutes). 
The time-development of the magnetic field after this time will 
be the subject of a planned future paper. 
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Fig. 1. Example of radial cut in NaDi filtergrams near the limb. 
The image on left shows the full field of view with the mag- 
netic feature of interest circled. On top right the smaller image 
is a zoom in of the feature and the plot below shows the radial 
cut. The center- side and limb- side amplitudes, Gc and ai are in- 
dicated. 



We make two sets of simulations with two diff'erent upper 
boundary conditions (BCs) for the magnetic field: vertical and 
potential field. In the potential filed case we considered both a 
case where the field was vertical at infinity and where it was 
slightly inclined. The inclination was found to make a negligible 
diff'erence for the expansion. Otherwise the BCs for the two sim- 
ulations are the same: velocity EC is open at the top and bottom, 
and all BCs are periodic in the horizontal directions. As for the 
thin flux tube slabs we produce synthetic observables at vary- 
ing LOS angles from snapshots of the simulation taken after the 
main transient phase due to the initial conditions is over (approx- 
imately 12 minutes in solar time). Since the simulations extend 
only to approximately 700 km above the t=1 surface, the ob- 
servables are formed in a region strongly influenced by the BCs 
(this is discussed in section 4.2). Rays entering at an inclined 
LOS are even more strongly influenced by this. These caveats as 
well as the assumption of LTE for the radiative transfer should 
be kept in mind when comparing with observations. 



3. Hinode observations 

3.1. NaDi filtergrams 

Due to expansion of magnetic flux with height the appearance of 
magnetic features in NaDi Stokes V filtergrams changes from 
unipolar at disk center to bipolar close to the limb (Fig. [T]). 
Additionally, due to foreshortening eff'ects the features become 
more elongated parallel to the limb. The size of the apparent 
bipolar features varies in the direction parallel to the limb from 
a few arcseconds to well over ten. In the radial direction the size 
(without correcting for foreshortening) is around four arcsec- 
onds at yu < 0.4. No significant further change in size is seen to- 
wards the limb. As one moves away from the disk center towards 
the limb the circular polarization radial cuts change from having 
a single peak (unipolar cuts) to having two peaks with opposite 
signs (bipolar cuts). The limb- side peak is usually smaller. The 
ratio of the peaks (j^^ see Fig.[T]) decreases toward the limb 
and will be used as a proxy for expansion. 
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Fig. 2. Histogram of number of uni- (dashed line) and bipolar 
(dotted line) circular polarization radial cuts in the NFI data as a 
function of yu. 



Histograms of the number of unipolar and bipolar radial cuts 
as a function of yu (Fig.[2l) display the change from unipolar to 
bipolar: at disk center nearly all cuts are unipolar with either one 
or two peaks. The double-peaked cuts are from features with 
complex, i.e., not circular or compact, shapes. At the limb all 
cuts are double-peaked with oppositely signed peaks, i.e., bipo- 
lar. The apparent bipolar features become dominant at yu < 0.5. 
Because of the insufliicient sample it is not possible to say what 
the statistics are in yu-range 0.5-0.7 (very few observations at 
those yu- values). 

In most apparent bipolar features (102 of 1 15) the center- side 
peak has a larger amplitude than the limb-side peak (Fig. [3]). The 
ratio of the center- to the limb- side peak amplitude strongly de- 
creases towards the limb. Some scatter in Fig. [3] could be due 
to the network magnetic fields being inclined in the direction 
of the LOS or away from it. However, observations of photo- 
spheric network magnetic fie lds indicate that the field inclination 
is sm all (less than 10 deg, Sanchez Almeida & Martinez Pilletl 
Il994i) making inclined fields an unlikely explanation for most of 
the observed scatter. No diff'erence in the ratio is seen for fea- 
tures at the polar limb and at the equatorial limb, i.e., no sig- 
nificant diff'erence is seen when using the ratio as a proxy for 
expansion in the mostly unipolar polar regions and the likely 
more heterogeneous equatorial limb. This suggests that the ex- 
pansion studied here is driven locally by the properties of the 
individual magnetic features, independent of the magnetic sur- 
roundings. The diff'erence between quiet Sun and coronal holes 
is most pronounced in the upper layers of the atmosphere. Since 
the observables are formed low in the atmosphere, it is not sur- 
prising that no diflference is seen with the current expansion di- 
agnostic. Wiegelmann & Solanki (2004) showed using magnetic 
field extrapolations that in the quiet Sun both small and large 
loops are present while in coronal holes the large loops are re- 
placed by open field lines. The far more numerous small loops 
are more likely to have an eff'ect on the expansion in the photo- 
sphere than the large loops. 

3.2. SP scans 

Bipolar features are also seen in SP circular polarization maps 
close to the limb (Fig.©. However, unlike the NaDi filtergrams, 
in the purely photospheric iron lines a significant number of 



Pietarila et al: Expansion of magnetic flux concentrations 



10,0 r 



o 

o 



o 
en 



KO- 



0.1 






0.0 



0.2 



0.4 



0.6 



0.8 



Fig. 3. Scatter plot of NaDi ratios as a function of yu. Thick solid 
line shows the mean ratio in 0.1 yu bins and the thin lines are the 
mean plus/minus standard deviation. 




Fig. 4. As Fig.[T]but for the averaged circular polarization signal 
in the 630 nm Fe lines. 



single-peaked cuts are found even at small yu- values (Fig. O. 
Furthermore, the flux concentrations close to the limb often look 
more diflfuse in Fe than in NaDi making it difl&cult to judge 
where one feature ends and another one begins. Apparent bipolar 
features outnumber unipolar features at yu < 0.3. 

As in NaDi the ratio of the center- side to the limb-side peak 
amplitude decreases towards the limb, although not so clearly as 
for NaDi . Also, the ratio in Fe is roughly a factor two larger than 
in NaDi. A direct comparison between the two is not strictly 
valid since the parameters used for measuring the circular po- 
larization are not the same (filtergram data with a specific fil- 
ter transmission profile in NaDi and integrated area of circular 
polarization over the Stokes V spectral profile for the Fe lines) 
nor do the measurements have the same signal-to-noise ratio. 
Additionally, the line formation characteristics, such as Lande 
g-f actors, are diflTerent. The Fe lines likely also sample weaker 
magnetic fields not visible in the NaDi NFI observables. 

The peak amplitude of unipolar features' radial cuts decrease 
towards the limb although the scatter is very large. At disk cen- 
ter (yu > 0.8) in 48 of 63 cases the circular polarization signal 
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Fig. 5. Histogram of number of uni- (dashed line) and bipolar 
(dotted line) cuts as a function of yu for the SP data. 
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Fig. 6. Scatter plot of the SP data ratios as a function of yu. Thick 
solid lines show the mean ratio in 0. 1 yu bins and the thin lines are 
the mean plus/minus standard deviation. The dash-dotted line is 
the mean (in 0.1 yu bins) of the NaDi ratio. 



is larger in 630.1 than 630.2 indicating the 630.2 line may be 
Zeeman saturated. For oflT-disk center cuts (yu < 0.8) only 3 out 
of 66 cases have a larger peak amplitude in 630. 1 than 630.2. The 
decrease in peak amplitude of unipolar features towards the limb 
is due to the line-of sight becoming less aligned with the axis 
of the magnetic flux tubes (assuming that the magnetic field in 
the larger magnetic flux concentrations sampled in this study is 
mainly vertical in the photosphere) and the increasing formation 
height combined with the decrease in field strength with height 
(Sol anki et al.lll99 6). The peak amplitude of unipolar cuts as a 
function of yu is similar to the center-side amplitude in the bipo- 
lar cuts. The total polarization (circular and linear) is larger in 
the bipolar features consistent with the unipolar features being 
smaller flux concentrations which do not expand as strongly as 
the bipolar features. Since we do not see as many unipolar fea- 
tures in NaDi (in fact none were included in the analysis) this 
indicates that either the unipolar features further expand with 
height (assuming the NaDi signal is formed higher than the Fe 
lines and appear bipolar in NaDi, or they are not strong enough 
to reach or to be visible at the formation height of the NaDi 
signal in the NFI data. Due to the larger width and lower Lande 
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Fig. 7. Radius of thin flux tubes as a function of height. Colors 
correspond from darkest to Hghtest to 140, 200, 400 and 600 km 
ro. Thick Hne is for 5o=1300 G and thin Hne for 1400 G. Shown 
are also facular (soHd Hne) and plage (dash-dotted line) models 
as internal atmospheres, but the diflTerence between the two is 
barely visible. 



^-factor of the NaDi line, together with the NFI not being as sen- 
sitive as SP, the NFI data may sample only the stronger, larger 
magnetic features. Note that the NaDi signal does not nece ssar- 
ily always come from a greater height (Ide Wijn et al.ll2009l) . 



4. Synthetic data 

4.1. Thin flux tube and sheet approximations 

Besides the geometry (tube or sheet) the choice of flux tube ra- 
dius at z=0 km, ro, is the most important factor determining the 
expansion of the thin flux tubes (Fig. [7]). Note that the values 
of ro considered here (140, 200, 400 and 600 km. z = cor- 
responds to T = 1 in the external atmosphere) are comparable 
to or larger than the pressure scale height, i.e., outside the strict 
validity range of the thin flux tube approximation. The absolute 
expansion is significantly stronger for ro=400 and 600 km than 
for the tubes with smaller radii (although the relative expansion, 
r(z)/r(0), is the same for all). Compared to the choice of radius, 
the choice of internal model (facular or plage) or ^o (magnetic 
field strength at z=0 km, 1300 or 1400 G) does not produce 
dramatically diff'erent expansion rates as shown in Fig. [TJ The 
plage and network models used here describe the network and 
plage regions as a whole, i.e., they do not specifically address 
flux tubes. The discre pancy between these mod els and flux tube 
models, such as, e.g. jBriand & Solankil (Il995h . is largest in the 
low photosphere where the temperature rise is stronger in the 
flux tube models. This leads to slightly stronger expansion of 
the tube in the lower layers. In the following, we show synthetic 
observables from slabs with the facular model as the internal at- 
mosphere and ^0=1300 G, where ^o corresponds to B in the the 
flux tube at = 0. Ratios resulting from diff'erent models (atmo- 
spheric model, ^o) are not significantly diff'erent from the ones 
shown here. 

The radial cuts in NaDi and Fe become more asymmetric to- 
wards the limb (Fig. (8]). Local maxima in the cuts are produced 
close to the boundaries of the flux tube due to the hot- wall eflfect. 
When the tube radius, ro, is comparable to the pressure scale 
height (140 km or 200 km) the tube does not expand enough or 
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Fig. 8. Synthetic radial cuts (Hinode resolution) from thin flux 
tubes at diff'erent yu, decreasing from top to bottom (ji is given 
within each of the frames in the left column) The various curves 
refer to diff'erent ro: 200 km in turquoise, 400 km in blue, and 
600 km in pink. In left column the averaged Fe 620 nm and in 
the right column Na Di Stokes V amplitudes (sign given by sign 
of blue lobe) are plotted. 



at heights low enough to produce clearly bipolar features in the 
Fe and NaDi observables as shown in Fig.[9l where the results 
are summarized for all tubes with ratios below ten. Only when 
the tube radius is > 400 km do bipolar features begin to appear at 
yu=0.5. The smallest ratios are from the 7?o=600 km case where 
the expansion of the magnetic field is the strongest. Using the 
amplitude instead of area results in slightly diflferent ratios. The 
diff'erences between area and amplitude are largest for the small- 
est yu-values, but both observables show the same trends. The 
eflfect of convolving and rebinning the observables to Hinode 
resolution only slightly alters the ratios. Normalizing the Stokes 
V profiles to the local continuum intensity prior to computing 
the observables (area, amplitude) slightly changes the ratios. The 
NaDi ratios are slightly higher than the Fe ratios, contrary to the 
observations. For all the observables (original resolution, Hinode 
resolution, area, amplitude, normalized or not) the trend is the 
same: tubes with ro=400 or 600 km produce bipolar features at 
yu-values below 0.5 and the ratios decrease towards the limb. The 
values of the resulting ratios are consistent with the observed 
ones. 

The radial extent of the apparent bipolar features in the flux 
tubes is at most ^3 arcseconds at disk center. The sizes are even 
smaller close to the limb. No large diff'erences are seen in width 
when comparing the synthetic Fe and NaDi cuts, but in general 
the features appear wider in Stokes V area cuts than amplitude 
cuts. This is due to the profiles near the flux tube boundaries 
having asymmetric, multi-lobed shapes. The integrated absolute 
area in these profiles can be relatively large whereas the ampli- 
tude tends to be fairly small. If the foreshortening eff'ect close 
to the limb is taken into account the synthetic features are sig- 
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Fig. 9. Amplitude (left) and area (right) ratios from synthetic ra- 
dial cuts produced from the thin flux tube models. Shown are 
only ratios smaller than ten. The amplitudes and areas are com- 
puted from Stokes V profiles which have not been normalized to 
the continuum intensity. 
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Fig. 10. As Fig.[7]but for the thin flux sheet models. 



nificantly smaller than those observed. In the disk center NFI 
data, diameters of network flux concentrations range on average 
from 1 to 3 Mm, i.e., the ro=600 km tube is still relatively small. 
Such small features as these resulting from the flux tube models 
would be very difl&cult to identify close to the limb. An impor- 
tant factor not included in the present analysis is that we expect 
the network to be composed of many smaller flux tubes/sheets, 
which may not expand like a single sheet/tube. We address this 
with the MHD simulations (next section). 

Fig. [TOl shows the radius of the thin flux sheets as a func- 
tion of height. As expected, the fanning out is stronger and takes 
place at lower heights in the sheets than in the tubes. Unlike 
the tubes, sheets with ro comparable to the pressure scale height 
show significant expansion with height. The resulting observ- 
ables (Fig.[TT]) diff'er from the tube case in that 1) the radial cuts 
become bipolar in appearance at larger yu-values 2) sheets with 
smaller ro produce bipolar cuts. Also the radial extent of the ap- 
parent bipolar features is larger in the sheets. The ratios (Fig. [12]) 
are similar to the tube case. The largest sheets produce bipolar 
features already atyu=0.77. The switch over from uni- to bipolar 
in the SOT observations takes place between /d=0.1 and 0.5. As 
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Fig. 11. As Fig.[8]but for the thin flux sheet models. 



in the tube case the NaDi ratios are slightly larger. Contribution 
and response functions for the synthetic profiles show that the 
Stokes V signals of the three lines are formed at a very similar 
height: the NaDi core is formed higher than the Fe line cores 
but the Stokes V signal is mostly from the line flanks which are, 
at least in LTE, mostly formed in the same height range. The 
only exception are rays entering the flux tube at an inclined an- 
gle: compared to the Stokes V signal of the Fe lines, the bulk 
of the NaDi Stokes V signal from the center-side of the tube 
is formed slightly higher, where also the line of sight magnetic 
field is slightly higher. An additional factor which may play a 
minor role is saturation. Since the Fe lines are more Zeeman 
sensitive they become Zeeman saturated at lower field strengths 
than NaDi. Zeeman saturation of profiles from the disk center 
side of the flux tube would lead to a smaller ratio. We did not 
explore this further nor did we test varying the microturbulent 
velocity in the models. Increasing it would broaden the lines and 
result in a reduced Stokes V amplitude. 



4.2. MuRAM simulations 

The magnetic field is highly inhomogeneous in the MuRAM 
simulations. In the vertical field BC simulation (Fig. [13]) the 
magnetic field remains confined to the strip in the middle of the 
simulation domain also after the transient phase. Even high up 
in the simulation domain the area around the magnetic strip re- 
mains void of magnetic field. The main magnetic structure is 
composed of distinct individual concentrations rooted in the in- 
tergranular lanes. These structures do not entirely merge at any 
height, but their expansion leads to a decrease in field strength 
with height and only horizontal fields in the individual con- 
centrations have a bipolar-like structure: a horizontal black and 
white pattern in the horizontal field, i.e., oppositely directed, in 
Fig.[T3]is visible only in the individual magnetic flux concentra- 
tions and not on a larger spatial scale. 
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Fig. 12. As Fig.[9]but for the thin flux sheet models. 

In contrast, in the potential field BC snapshot (Fig. \T^ the 
magnetic field fills the entire upper part of the domain, forming 
a large scale canopy-like configuration. Individual flux concen- 
trations are still present close to t= 1 , but towards the top they 
merge into the canopy- structure. The canopy structure is very 
inhomogeneous on small scales and the eflTect of down flows re- 
lated to convective overshooting is clearly seen. Portions of the 
canopy are dragged down to the bottom of the photosphere by 
the flows: there is a significant amount of magnetic field, both 
vertical and horizontal, in the intergranular lanes at z=0 km. A 
more detailed view of the potential field BC simulation is given 
in Fig {151 A large portion of the field outside the strip is horizon- 
tal even at low heights, but it is not as homogeneous as higher 
up. The relative (to the vertical component) strength of the hor- 
izontal field increases with height as is seen in the cross-cuts of 
the average total and horizontal magnetic field strengths. This is 
caused by both the vertical component becoming weaker with 
height but also by the horizontal component becoming stronger. 

The synthetic observables shown in Figs. [161 and [T7] give 
an impression of how observables from such magnetic config- 
urations would look like. We have not added ordinate scales 
to the synthetic radial cut figures because the main diagnostic 
(i.e., the ratio of the diskward to the limb ward peak amplitude of 
Stokes V) depends on the expansion and disk position, but not 
on the amplitude of V. Since the magnetic field is inhomoge- 
neous at all heights, the resulting cuts exhibit small-scale struc- 
ture. Since the spectral lines are formed at a roughly constant 
optical depth instead of geometric height, the strong changes in 
density across the simulation domain (flux concentrations and 
the surroundings) lead to a strongly varying geometric forma- 
tion height range. The radial cuts do not mimic the cross-cuts of 
the magnetic field shown in Fig. \T5\ After taking into account 
the spectral and spatial resolution of Hinode the cuts from the 
vertical BC simulation do not at any LOS viewing angle appear 
bipolar. In contrast, cuts from the potential BC simulation dis- 
play apparent bipolar features at ji >0.34. This is closer to the 
limb than in the observations or thin flux sheet and tube models 
(yu = 0.5). The ratio decreases with decreasing//- value. However, 
the ratios remain fairly large, ^5. 

5. Discussion 

The SOT observations show clearly the expansion of magnetic 
flux concentrations with height as a change from unipolar radial 
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Fig. 13. Snapshots of the magnetic field in the simulation with 
a vertical field at the top boundary. Top: horizontal cut through 
the simulation domain showing the x-componentof the magnetic 
field at a geometric height of roughly 400 km above the r = 1 
level. Dotted line shows the location of the vertical cuts shown in 
the lower frames. Middle and bottom: vertical cuts of the x- and 
z-component of the magnetic field. Dotted line shows the height 
of the horizontal cut shown in the top panel. 



cuts at disk center to bipolar near the limb. Using the ratio of the 
radial cuts' center- to the limb-side peak as a proxy for expan- 
sion indicates that the magnetic field has fanned out more at the 
formation height of the NaDi NFI signal than at the formation 
height of the Fe I 630 nm SP signal. The NaDi ratios are smaller 
than the Fe I and the switch over from unipolar to bipolar cuts 
takes place closer to disk center, i.e., at yu=0.5 in NaDi rather 
than at ji ^0.3 in Fe I. Furthermore, unlike the SP data the NFI 
data exhibit exhibit only very few unipolar features close to the 
solar limb. The apparent expansion between the SP and NFI data 
and the lack of unipolar features close to the limb in NFI data can 
be caused by several factors: the NaDi signal is formed higher 
and the features that would appear unipolar in SP look bipolar 
in NFI due to increased expansion of the field with height, as 
expected from flux tube models. Alternatively, the unipolar fea- 
tures are not strong enough or do not reach high enough to be 
visible in the NFI data. The line formation of the NaDi and Fe 
I 630 nm lines are somewhat diflferent, e.g., the eff'ective g value 
of the NaDi line is smaller than for either Fe line making it less 
sensitive to low magnetic field strengths or fluxes. The diff'erence 
may be accentuated by other factors such as the NFI data being 
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Fig. 14. As Fig.[T3]but for the simulation with a potential field at 
the top boundary. 



filtergram while the Fe I signals are deduced from spectrograph 
data. The diff'erent noise levels and pixel sizes of the instruments 
may also play a role. To fully compare the two one should use 
simultaneous NFI and SP observations at various disk positions. 

Expansion based on the observed ratios and ratios derived 
from the thin flux sheet and tube approximations are in good 
agreement. The synthetic ratios are of similar magnitude and the 
switch over from unipolar to bipolar takes place at similar [i- 
values, i.e., around 0.5. Since the sheets expand faster, the thin 
sheet model results in bipolar features closer to disk center, first 
measurable ratios appear atyu ^0.8, and at smaller ro: 200 km for 
the sheet and 400 km for the tube models (assuming the observa- 
tions are perpendicular to the elongated direction of the sheet). 

Unlike in the observations, in the thin sheet/tube models the 
ratios between the limb-to-center side Stokes V amplitudes and 
areas of the Fe lines are smaller than the NaDi ratios. While the 
NaDi line core is formed higher than the Fe line cores, the line 
flanks, from which the bulk of the Stokes V signal is from, are 
formed at roughly the same height in the sheet and tube models 
(based on response functions computed in LTE). The diflTerence 
between the synthetic and observed ratios is probably at least 
partly due to the radiative transfer being done in LTE and a sam- 
pling bias due to the diflTerences in the SP and NFI data. Also, 
factors such as model magnetic field strength and microturbu- 
lence may influence the ratios. If the disk center-side of the flux 
concentration has a strong enough field it may become Zeeman 
saturated before the limb-side, leading to a relatively reduced 
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Fig. 15. Simulation with a potential field at the top boundary. Top 
left panel: average (over y-direction) x-component of magnetic 
field (top) and average magnetic field strength (2 lower panels) 
at three diff'erent heights. Top right: histograms of B^ at three 
diff'erent heights. Color code for top panels is given by the leg- 
end. Middle row: horizontal cuts of magnetic field inclination 
at three diflferent heights. Shown is l/3rd of the domain size for 
each height. Bottom: as middle row but for the magnetic field 
azimuth. 



center side amplitude and a decreased ratio. Since the Fe lines 
have a larger Zeeman sensitivity this eff'ect would first take place 
in the Fe lines. 

The observations clearly show the expansion of the field 
which is found to be similar to the expansion derived from us- 
ing the thin flux sheet and tube approximations. This result is in 
agreement with lSolanki et al.l (Il999l) . In other aspects the mag- 
netic network in the Sun is probably not well described by a 
simple flux tube or sheet scenario. In reality the network is dy- 
namic and likely composed of individual flux concentrations 
which merge at some height. Observati ons of chromospheric 
lines, such as H-or (e.g.Jde Pontieu et aP 2007), the Ca K (e.g., 
IZirinlll974 iPietarila et al.ll2009h and Ca II infrared triplet lines 
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Fig. 16. Synthetic radial cuts from the vertical field BC simula- 
tion at diff'erent yu- values (given in the left frames). Left Fe and 
right NaDi amplitudes (profiles are not normalized to the con- 
tinuum intensity). Gray line is at original simulation resolution 
and black at Hinode resolution. The LOS angle is indicated by 
the arrow in each figure (on the left). 
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Fig. 17. As Fig.[T6lbut for the potential field BC simulation. 



around 850 nm (e.g.. IVecchio et al.ll2007l) . exhibit fibril struc- 
tures which probably outline magnetic field lines that are more 
heated than others, i.e., the fibrillar canopy is thermally inhomo- 
geneous. The magnetic field in contrast is likely quite smooth as 
indicated by the homogeneous appearance of both the uni- and 
bipolar patches in the SOT data. 

The choice of upper BC for the magnetic field in the MHD- 
simulations has a significant eff'ect on the expansion properties of 
the field. The BC aff'ects not only the upper portion of the simula- 
tion domain but eff'ects are also seen at the t= 1 level, namely the 
existence of horizontal and vertical field also outside the initial 
magnetic strip. In the potential field BC simulations horizontal 
magnetic field is present everywhere and a canopy-like structure 
is formed at the top of the domain whereas in the vertical field 
BC simulation the field remains confined in the strip. We plan to 
address in a later paper the time evolution of the magnetic flux in 
the potential field BC simulation to study in detail the transport 
of magnetic field from the top of the domain, i.e., the canopy, to 
the bottom of the photosphere. Note that with the term canopy 
we are considering the canopy-like expansion of the field, not the 
actual chromospheric canopy. The choice of BC aflTects strongly 
the synthetic radial cuts: no apparent bipolar features are seen in 
radial cuts from the vertical BC simulation. In contrast, a bipolar 
feature begins to appear at jj = 0.34 in the potential BC case. 
The switch over from unipolar to bipolar takes place closer to 
the limb and the resulting ratios are larger than those seen in the 
observations or in the zeroth order thin flux tube/sheet models, 
indicating that the field does not fan out enough. 

A stronger expansion in the potential field BC simulation 
could be achieved by, e.g., having two strips of opposite mag- 
netic polarities in the initial setup or increasing the simulation 
domain size horizontally. The former would result in a canopy 
structure connecting the opposite polarities. However, based on 
the NFI ratios, the magnetic field in the mostly unipolar polar 
regions expands at a similar rate as in the equatorial limb re- 
gions, where mixed polarity fields are found. The latter alter- 
native, larger simulation domain, would result in stronger fan- 
ning out because of the periodic BCs in the horizontal direc- 
tion. Currently the domain size is 24 Mm which is comparable 
to the typical size of internetwork cells (30 Mm, B eckerslll968h . 
Therefore a larger domain would not necessarily be more re- 
alistic. Placing the upper boundary higher might also result in 
stronger expansion since the field would have more volume to 
expand in. However, the current treatment of radiative transfer 
in LTE would become increasingly incorrect. Finally, the field 
may not yet have reached a potential state at the upper boundary, 
so that even this boundary condition may not be appropriate. 

We plan to extend the current study by using spectropolari- 
metric observations of photospheric and chromospheric lines to 
study in detail the network magnetic fields. 

5.1. Conclusions 

The circular polarization signal of magnetic flux concentrations 
is known to change from unipolar at disk center to bipolar 
near the limb, consistent with the magnetic field fanning out 
with height. Solanki et al. (1999) showed with observations then 
available that the relative expansion rates of features of vari- 
ous sizes, were consistent with rates derived for simple thin flux 
tubes. The data from the Hinode satellite, with their high spa- 
tial resolution, low noise and freedom from atmospheric seeing 
eflfects, provide a new opportunity to address the question with 
better statistics and with an emphasis on the variation from disk 
center to the limb. We again found the observations to be in good 
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agreement with expansion properties derived from the thin flux 
sheet and tube approximations, confirming the previous results 
with the new data and a more extensive analysis. Since the sig- 
nal we are interested in forms higher in the atmosphere near the 
limb, we found that realistic numerical simulations, which ex- 
tend only to 700 km above logirsoonm) = 0, are less useful in 
modeling these observations, with the results strongly depend- 
ing on the choice of boundary conditions for the magnetic field. 
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